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BESCHREIBUNG 

Verfahren zur elastographischen Untersuchung eines Objektes 

Die Erfindung bezieht sich auf ein Verfahren zur elastographischen Untersuchung eines 
Objektes sowie auf eine Anordnung zur Durchfiihrung eines solchen Verfahrens und auf 
5 ein Computerprogramm zur Steuerung einer solchen Anordnung. 

Bei der aus der US-Appln. 09/743659 (PHD 99-152) bekannten MR-Elastographie (MR 
= Magnetresonanz) wird die Tatsache ausgenutzt, dass die Phase in einem MR-Bild des 
Objektes sich in Folge von mechanischen Schwingungen andert, die auf das Objekt 

10 einwirken. Das Ausmafl dieser Anderung hangt von der durch die mechanischen 
Schwingungen verursachten Auslenkung des Gewebes ab, Somit lasst sich aus MR- 
Phasenbildern, d. h. Bildern, die die Phase der Kernmagnetisierung darstellen, 
Information uber bestimmte mechanische Parameter des Gewebes ableiten, z. B. uber die 
Elastizitat Bei der Ultraschall-Elastographie ist sogar eine direkte Messung der 

15 Auslenkung moglich. 

Bosartige Lasionen wie z. B. ein Mamakarzinom unterscheiden sich durch ihre Elastizitat 
deutlich von normalem gesunden Gewebe. Doch gibt es auch gutartige Lasionen, die 
ahnliche elastische Eigenschaften haben. Anomale elastische Eigenschaften sind daher fur 
20 sich alleine noch kein hinreichendes Indiz fiir die Malignitat eines Gewebes. 

Aufgabe der vorliegenden Erfindung ist es, ein elastographisches Verfahren so 
auszugestalten, dass sich zusatzliche Informationen ergeben, die eine bessere 
Unterscheidung zwischen unterschiedlichen Gewebetypen ermoglichen. Diese Aufgabe 
25 wird erfindungsgemaB gel6st durch ein Verfahren zur Untersuchung eines Objektes mit 
den Schritten; 
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a) Anregung des Objektes mit mechanischen Schwingungen, 

b) Messung des zeitlichen Verlaufs der Auslenkung von Voxeln des Objektes, die den 
von den Schwingungen im Objekt verursachten Wellen ausgesetzt sind, aus ihrer 
Ruhelage, 

5 c) Bestimmung der nichtlinearen Verzerrungen aus dem zeitlichen Verlaufs der 
Auslenkung, 

d) Auswertung der nichtlinearen Verzerrungen. 

Die Erfindung basiert auf folgenden Uberlegungen und Erkenntnissen: 
10 Wenn der Druck auf biologisches Gewebe verandert wird, andert sich die Deformation 
des Gewebes nicht in gleichem MaBe wie der Druck, weil der Elastizitatsmodul zunimmt 
(das Gewebe wird harter). Dieser nichtlineare Zusammenhang zwischen Druck und 
Deformation ist bei malignen Lasionen besonders ausgepragt und hat zur Folge, dass bei 
einer Anregung des Gewebes durch mechanischen Schwingungen der zeitliche Verlauf 
15 der Auslenkung nicht mehr dem zeitlichen Verlauf der mechanischen Schwingungen 
entspricht, sondern (nichtlineare) Verzerrungen aufweist Daher sind die nichtlinearen 
Verzerrungen im zeitlichen Verlauf der Auslenkung eines Voxels ein MaB fur die 
Nichtlinearitat, das sich fur die Klassifizierung Gewebes eignet. Diese Verzerrungen 
fuhren bei einer Anregung mit zeitlich beispielsweise sinusfbrmig verlaufenden 
20 mechanischen Schwingungen dazu, dass die Auslenkung der Voxel in dem Objekt nicht 
mehr rein sinusfbrmig erfolgt, sondern zusatzlich hohere Harmonische der mechanischen 
Schwingungen aufweist - insbesondere solche mit der doppelten Frequenz. 

Man kQnnte jedes Voxel einzeln erfassen und die dafur ermittelten nichtlinearen 
25 Verzerrungen (und gegebenenfalls noch die Elastizitat) ermitteln. Die Ausgestaltung 
nach Anspruch 2 erleichtert dem Untersucher aber die Diagnose fur einen groBeren 
Bereich. 



30 



Anspruch 3 gestattet die Durchfiihrung der Erfindung mit den Mitteln der MR- 
Elastographie, wodurch es moglich wird, neben der Nichtlinearitat auch cllen 
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Elastizitatsmodul ortsabhangig zu bestimmen. Grundsatzlich kann das Verfahren aber 
auch mit Hilfe von Ultraschall ausgefuhrt werden. 

Anspruch 4 beschreibt eine Ausgestaltung des Verfahrens nach Anspruch 3, die es 
5 ermoglicht, gleichzeitig Auslenkungen eines Voxels mit der Grundfrequenz und mit dem 
Doppelten der Grundfrequenz zu ermitteln. Stattdessen ware es aber auch moglich, das 
magnetische Gradientenfeld aus einer Komponente mit der Grundfrequenz und einer 
Komponente mit dem doppelten der Grundfrequenz zusammenzusetzen. 
Die weitere Ausgestaltung nach Anspruch 5 verhindert den StoreinfluB auf die Phasen 
10 aufgrund von Feldinhomogenitaten. In diesem Fall miissen die beiden einzelnen 
sinusformigen Schwingungen des Gradientenfeldes dieselbe Phasenlage haben. 

Die Ausgestaltung nach Anspruch 6 tragt der Tatsache Rechnung, dass mit der 
(ortsabhangigen) Nichtlinearitat und Elastizitat das Gewebe besser klassifiziert werden 
15 kann. Bei dem Verfahren nach Anspruch 3 sind hierfur aber keine zusatzlichen 
Messungen erforderlich. 

Anspruch 7 definiert ein MR-Gerat zur Durchfiihrung des erfindungsgemaBen 
Verfahrens, und Anspruch 8 beschreibt ein Computerprogramm zur Steuerung eines 
20 solchen MR-Gerates. 

Die Erfindung wird nachstehend anhand der Zeichnungen naher erlautert. Es zeigen: 
Fig. 1 ein Diagramm zur Einordnung unterschiedlicher Gewebetypen, 
Fig. 2 ein MR-Gerat zur Durchfiihrung des erfindungsgemSBen Verfahrens, 
25 Fig. 3 ein Ablaufdiagramm fur die Acquisition von MR-Phasenbildern, 
Fig. 4 eine der zur MR-Akquisition dienenden Sequenzen, 
Fig. 5 eine schematische Darstellung eines Tetis des Auswerteverfahrens und 
Fig. 6 ein Ablaufdiagramm des Auswerteverfahrens. 



30 
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Zwischen dem Druck P auf eine Probe und der unter diesem Druck erfolgenden 
Verformung e der Probe besteht bekanntlich die Beziehung 



P = sE 



(1) 



5 

Dabei ist E die Elastizitat des Gewebes. Fur die Verformung e gilt die Beziehung 



10 Dabei ist L 0 die Ausdehnung der nicht komprimierten Probe und L die Ausdehnung der 
komprimierten Probe. Wie eingangs erlautert, ist bei Proben aus biologischem Gewebe 
die Verformung s dem Druck P nicht proportional, weil die Elastizitat mit zunehmender 
Deformation zunimmt. Diese Zunahme kann in guter Naherung beschrieben werden 
durch die Gleichung 



Dabei ist Eo die Anfangselastizitat (d. h. die Elastizitat ohne Kompression), und der 
Exponent a ist ein Mafl fiir die Nichtlinearitat. FOr a = 0 ergibt sich ein linearer 
20 Zusammenhang zwischen P und e. Je grSBer a ist, desto starker nichtlinear wird dieser 
Zusammenhang. 

Fig. 1 zeigt ein Diagramm mit der Elastizitat Eo als Ordinate und dem die Nichlinearitat 
kennzeichnenden Exponenten a als Abszisse. Jeder Probe lasst sich in Fig. 1 ein Punkt 
25 zuordnen, der die Elastizitat E 0 und den Exponenten a fiir diese Probe kennzeichnet. In 
dem Diagramm sind eine Anzahl von Gewebeproben durch je ein Symbol markiert Mit 
einem „+" sind Gewebeproben aus Mammakarzinomen bezeichnet. Ein „o" bezeichnet 
gutartige Lasionen in Form eines Fibroadenoms, ein „□" bezeichnet Fettgewebe und 
durch ein „A" sind Gewebeproben von Wucherungen in den Milchgangen bezeichnet, 



s = (L 0 -L)/L 0 



(2) 



15 



E = E 0 e' 



oce 



(3) 
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die von einem gutartigen in einen bSsartigen Zustand ubergehen konnen. 

Gewebeproben von Mamakarzionomen haben eine groBere Elastizitat (d. h. sie sind 
harter) und gleichzeitig einen starker nichtUnearen Zusammenhang zwischen dem Druck 
5 P und der Deformation 8 als gesundes Gewebe. Man kann Gewebeproben daher nach 
diesen Parametern klassifizieren. Gewebeproben rechts von der in das Diagramm 
eingezeichneten Geraden deuten auf eine maligne Lasion hin, wahrend Gewebeproben 
links davon ein Zeichen flir eine gutartige Lasion sind. 

10 Die Einordnung in das Diagramm nach Fig. 1 basiert auf in-vitro-Untersuchungen von 
Gewebeproben, bei denen die Gewebeproben einem steigenden Druck P ausgesetzt 
wurden und die Verformung e gemessen wurde. Das Gewebe von lebenden Patienten 
lasst sich auf diese Weise nicht ohne weiteres, untersuchen. Die Erfindung beruht auf der 
Erkenntnis, dass man dieselben Parameter nicht invasiv an lebenden Patienten mit 

1 5 Mitteln der Elastographie bestimmen kann. Wenn ein lebendes Objekt durch stationare 
mechanische Schwingungen angeregt wird, dann musste die Auslenkung der Voxel in 
dem mechanischen Gewebe den gleichen zeitlichen Verlauf haben wie die mechanischen .. ; 
Schwingungen, wenn die Elastizitat konstant und von der Verformung unabhangig ware. ■ 
Befindet sich ein Voxel jedoch innerhalb eines Gewebes, das eine nichtUnearen 

20 Zusammenhang zwischen dem Druck P und der Verformung e aufweist, dann kann die 
Auslenkung des Voxels dem zeitlichen Verlauf der Schwingungen nicht mehr genau 
folgen und es entstehen Verzerrungen, die ein MaB fur die Nichtlinearitat sind. Die 
Messung dieser Verzerrungen bzw. Nichtlinearitat wird im Folgenden am Beispiel der 
MR-Elastographie erlautert. 

25 

Fig. 2 stellt schematisch einen Querschnitt durch ein dafiir geeignetes MR-Gerat dar 
sowie die zu dessen Betrieb erforderlichen Schaltungs- und Software-Komponenten. Das 
MR-Gerat umfasst einen Hauptfeldmagneten 30 zur Erzeugung eines horizontal und 
parallel zur Zeichenebene der Fig. 2 verlaufenden Hauptmagnetfeldes B 0 . Im Innern 
30 dieses zylinderfbrmigen Hauptfeldmagneten befindet sich ein Satz 3 1 von 
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Gradientenspulen zur Erzeugung von Magnetfeldern mit Gradienten in drei zueinander 
senkrechten Richtungen. Im Innern dieser Gradientenspulen wiederum befindet sich eine 
Hochfrequenzspule 32 zur Erzeugung eines Hochfrequenz-Magnetfeldes. 

5 Die Strome fur die Gradientenspulen werden von einer Generatoranordnung 33 geliefert, 
wahrend die StrSme fur die Hochfrequenzspule 32 von einem Hochfrequenzgenerator 34 
geliefert werden. Der zeitliche Verlauf der von den Generatoren 33 und 34 erzeugten 
StrSme wird von einer Steuereinheit 35 vorgegeben und gesteuert. Die bei der 
Untersuchung in den Mammae entstehenden MR-Signale werden von nicht naher 

10 dargestellten MR-Spulen empfangen und in einem Empfanger 36 weiterverarbeitet. Die 
Ausgangssignale des Empfangers 36 werden einer Bildverarbeitungseinheit 37 zugefuhrt, 
die aus den empfangenen MR-Signalen ein MR-Bild rekonstruiert, das auf einem 
Monitor 38 wiedergegeben wird. 

1 5 In dem durch die Spulen 3 1 und 32 defmierten Untersuchungsbereich befindet sich eine 
Patientin 47, die auf einem Kissen 41 gelagert ist, das auf einem Patiententisch 2 Uegt. 
Vor dem Kissen 41 ist ein Mammographiezusatz 1 angeordnet, der die gleiche H6he hat 
wie das Kissen 41 . Die Mammae der Patientin hangen in Ofihungen des 
Mammographiezusatzes hinein, wo sie unter leichter Kompression mechanischen 

20 Schwingungen ausgesetzt sind, die in horizontaler Richtung (vom FuB zum Kopf) und 
zeitlich sinusffcrmig verlaufen. Der zeitliche Verlauf der in dem Mammographiezusatz 1 
erzeugten mechanischen Schwingungen wird von der Steuereinheit 35 mit dem zeitlichen 
Verlauf der von den Generatoren 33 und 34 erzeugten Strdme synchronisiert. 

25 In der deutschen Patentanmeldung 10156178.4 (PHDE 010302) ist dieses MR-Gerat im 
Einzetnen beschrieben, worauf zur Vermeidung von Wiederholungen verwiesen wird. 

Fig. 3 erlautert den zeitlichen Ablauf bei der Akquisition der MR-Daten fur eine MR- 
Elastographie-Untersuchung mit dem MR-Gerat nach Fig. 2. Nach der Initialisierung 
30 (100) und nach der Lagerung der Patientin 40 entsprechend Fig. 1 werden deren 
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Mammae leicht komprimiert, so dass sie fixiert und weiterhin elastisch verformbar sind. 
Danach werden in dem Mammographiezusatz 1 zeitlich sinusformige mechanische 
Schwingungen erzeugt, die auf die Mammae der Patientin einwirken und die wahrend 
der gesamten MR-Akquisition anhalten (Block 101). 

Nachdem diese Schwingungen stationar geworden sind, wird zunachst im Block 102 
eine erste MR-Sequenz erzeugt, deren zeitlicher Verlauf in Fig. 4 dargestellt ist. Die 
Sequenz umfasst einen 90°-Hochfrequenzimpuls RFi und danach einen refokussierenden 
180°-Hochfrequenzimpuls RF 2 (4. Zeile von Fig.4). Zwischen den beiden 
Hochfrequenzimpulsen RF, und RF 2 wird eine einzelne, zeitlich sinusformig verlaufende 
Schwingung Gi des magnetischen Gradientenfeldes erzeugt (3. Zeile von Fig.4), deren 
Periodendauer genau halb so grofi ist (z. B. 10 ms bei einer Frequenz der mechanischen 
Schwingung von 50 Hz) wie die Periodendauer der mechanischen Schwingungen. Nach 
dem zweiten Hochfrequenzimpuls wird eine zweite Schwingung G 2 des magnetischen 
Gradientenfeldes G erzeugt, die exakt d en gleichen zeitlichen Verlauf hat wie die 
Schwingung G x und deren Beginn gegenuber dem Ende der ersten Schwingung urn eine, 
viertel Periode (5 ms) der mechanischen Schwingung versetzt ist. 

Wahrend der beiden Hochfrequenzimpulse RFi und RF 2 ist ein in z-Richtung 
verlaufendes magnetisches Gradientenfeld wirksam, das die Hochfrequenzimpulse nur in 
einer zur Langsrichtung der Patientin senkrechten Schicht wirksam werden lasst und das 
in Fig. 4 ebensowenig dargestellt ist wie der Phasenkodiergradient und der Lesegradient. 

In der ersten Zeile von Fig. 4 ist ein Ausschnitt aus dem zeitlichen Verlau f der 
mechanischen Schwingungen dargestellt und in der zweiten Zeile eine hohere 
Harmonische dieser Schwingungen mit der doppelten Frequenz. Solche hoheren 
Harmonischen kOnnen sich fortlaufend in Bereichen ergeben, in denen ein nichtlinearer 
Zusammenhang zwischen Deformation und Elastizitat besteht, wenn die durch die 
mechanischen Schwingungen angeregten Wellen sich in der Mama ausbreiten und auf 
diese Bereiche trifft. AuBerdem ist in den beiden ersten Zeilen durch die punktierten 
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Linien angedeutet, dass die mit der jeweiligen Frequenz ausgelenkten Kernspins durch 
den Impuls RF2 refokusiert sind. 

Durch die Gradientenschwingungen Gi und G 2 wird die Phase der durch die 
5 mechanischen Schwingungen ausgelenkten Kernspins geandert. Das AusmaB der 

Anderung hangt von der Starke der Auslenkung in Richtung des jeweiligen Gradienten 
ab, wahrend das Vorzeichen der Anderung von der Phasenlage zwischen der 
Gradientenschwingung und der Auslenkung der Kernspins abhangt. 

10 Man erkennt anhand der Fig. 4, dass das Vorzeichen der Gradientenschwingung Gi 
wechselt, wahrend das Vorzeichen der mechanische 50 Hz-Schwingung (1. Zeile) 
unverandert bleibt. Es ergibt sich daher keine Phasenanderung. Hingegen bewirkt die . 
zweite Gradientenschwingung G 2 eine Phasenanderung, weil die 50 Hz-Schwingung im 
selben Moment ihren Nulldurchgang hat wie die Gradientenschwingung G 2 . Bei den mit 

15 der Harmonischen (100 Hz) der Grundfrequenz ausgelenkten Kernspins bewirken beide 
Gradientenschwingungen Gi und G 2 gleich groBe Phasenanderungen, die sich aufgrund 
der dargestellten Verlaufe addieren. 

Die in Fig. 4 dargestellte Sequenz ist also sowohl fiir die Auslenkung mit der 
20 Grundfrequenz als auch fur die Auslenkung mit der doppelten Grundfrequenz sensitiv. 
Ware dies nicht der Fall, mussten die Grundfrequenzen der Auslenkung und deren 
Harmonische getrennt erfasst werden, wodurch die Messzeit verlangert wurde. 

Nachdem im Schritt 102 das resultierende MR-Signal empfangen worden ist, wird im 
25 Schritt 103 die MR-Sequenz fur eine andere, in z-Richtung (die in Fig. 2 horizontal 
verlauft) versetzte Schicht wiederholt (Block 103). Dabei wird lediglich die 
Tragerfrequenz der Hochfrequenzimpulse RFi und RF 2 so verandert, dass die andere 
Schicht angeregt wird, wahrend die zeitlichen Verlaufe der magnetischen 
Gradientenfelder unverandert bleiben. Die Schleife mit den Schritten 102 und 103 wird 
30 dann so oft durchlaufen, wie Schichten in dem zu untersuchenden dreidimensionalen 
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Bereich vorhanden sind, z. B. 20 mal, wobei jedesmal eine andere angeregt wird. 

Danach wird der Phasenkodierungsgradient im Schritt 104 geandert, und alle Schichten 
werden erneut angeregt und die darin erzeugten Spinechosignale erfasst. Jede einzelne 

5 Schicht wird dabei mit derselben zeitlichen Lage der Hochfrequenzimpulse und der 

magnetischen Gradientenfelder in Bezug auf die mechanische Schwingung angeregt. Der 
Phasenkodierungsgradient wird so oft geandert wie Phasenkodierungsschritte 
erforderlich sind, z. B. 128 oder 256 mal. Danach sind die MR-Daten von samtlichen 
Schichten akquiriert, so dass daraus ein MR-Bild rekonstruiert werden kann, das die 

10 Phase der Kernmagnetisierungsverteilung in dem durch die Schichten gebildeten 
dreidimensionalen Bereich darstellt. 

Nachdem die MR-Daten fur das erste MR-Phasenbild des dreidimensionalen Bereichs 
akquiriert worden sind, erfolgt im Schritt 105 eine zeitliche Verschiebung zwischen den 

15 MR-Sequenzen und den mechanischen Schwingungen, und zwar so, dass der Beginn 
eines MR-Experiments nunmehr urn einen Bruchteil einer Periode AT - z. B. AT = T/8 - 
versetzt ist (Schritt 105). Danach wird die Schleife mit den Schritten 102, 103 und 104 
erneut durchlaufen, so dass sich ein weiteres MR-Phasenbild von dem untersuchten , 
dreidimensionalen Bereich ergibt. Die beiden MR-Phasenbilder unterscheiden sich nur in 

20 den Bereichen, in denen das Gewebe bzw. die dort angeregten Kernspins in x-Richtung 
ausgelenkt worden sind. Danach werden weitere Phasenbilder erzeugt, wobei der 
zeitliche Versatz zwischen den mechanischen Schwingungen und der Sequenz verandert 
wird (z. B. AT = 2T/8, 3T/8 ... 7T/8). Die vom Benutzer eingegebenen oder durch seine 
Eingabe aufgerufenen Daten kdnnen auf dem Monitor 8 oder aber einem gesonderten 

25 Display wiedergegeben werden. 

SchlieBlich sind die MR-Daten von einer Anzahl von MR-Phasenbildern (im Beispiel 8) 
akquiriert, die von der Auslenkung der Kernspins in x-Richtung abhangen- wenn 
angenommen wird, dass der Gradient der Gradientenfelder Gi, G 2 in x-Richtung verlauft. 
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Urn auch die Auslenkung in einer anderen Richtung als der x-Richtung erfassen zu 
konnen, werden die sinusformigen Gradienten Gi und G 2 statt in x-Richtung in einer 
anderen Richtung erzeugt (Block 106), z. B. in y-Richtung, so dass - nachdem der 
beschriebene Durchlauf der Schleifen mit den Schritten 102, 103, 104 und 105 vollendet 

5 ist - ein Satz von MR-Phasenbildern des dreidimensionalen Bereichs zur Verfugung 
steht, der von der Auslenkung der Kernspins in y-Richtung beeinflusst ist. SchlieBlich 
wird dies fur z-Richtung wiederholt (d. h. die sinusfbrmigen Gradienten Gi und G 2 
werden in z-Richtung appliziert), so dass schlieBlich drei Satze von je 8 MR- 
Phasenbildern zur Verfugung stehen, die von der Auslenkung der Kernspins in x,- y, 

10 oder z-Richtung abhangig sind. Danach ist das Akquisitionsverfahren beendet (Block 
107). 

Anhand der Fig. 6 wird nunmehr die Verarbeitung der akquirierten MR-Daten erlautert. 
Nach der Initialisierung im Block 200 werden im Schritt 201 - wie schon erlautert - aus 
15 den akquirierten MR-Daten MR-Phasenbilder erzeugt. Im Schritt 202 wird fur jede der 
drei Raumrichtungen der zeitliche Verlauf der Auslenkung u x , u y und u z fur die einzelnen 
Voxel des dreidimensionalen Untersuchungsbereichs ermittelt. 

Zur Erlauterung wird auf die schematische Darstellung in Fig. 5 verwiesen. In der oberen 
20 Zeile sind die drei Satze von MR-Phasenbilder schematisch dargestellt, die von der 

Auslenkung in x-, y- bzw. z-Richtung abhangen. In jedem der zu einem Satz gehorenden 
MR-BUder wird die Phase in dem selben Voxel betrachtet. Aus den acht Bildern ergeben 
sich acht Stutzstellen fur den zeitlichen Verlauf der Phase bzw. der Auslenkung in x-, y- 
oder z-Richtung in dem betreffenden Voxel. Wenn dieses Voxel in einem Bereich mit 
25 nichtlinearen elastischen Eigenschaften liegt, verlauft die Auslenkung in diesem Voxel 
zeitlich nicht rein sinusformig, weil die Auslenkung in diesem Voxel neben einer 
Komponente mit der Frequenz der Grundschwingung auch noch hdhere Harmonische 
enthalt, insbesondere eine Komponente mit der doppelten Grundfrequenz. 
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Zur Ermittlung der Amplitude dieser Komponenten wird im Schritt 203 eine Fourier- 
Transformation uber die durch acht Stutzstellen definierte Phase bzw. Auslenkung des 
Voxels durchgeftthrt. Aus der Fourier-Transformation resultieren drei GroBen ai x , a Xy 
und ai Z) die die Amplituden der Auslenkung in x-, y- und z-Richtung mit der 
5 Grundfrequenz darstellen, und drei weitere GrfiBen ^ a 2y und aa* die die Amplituden 
der Auslenkung in x-, y- und z-Richtung mit der doppelten Grundfrequenz darstellen. 
Daraus lasst sich die Amplitude a, der Komponente mit der Grundfrequenz berechnen 
nach der Gleichung 



10 a, =«JaI+4+4 

Fiir die Amplitude der Komponente mit der doppelten Grundfrequenz gilt 



15 

Setzt man die beiden Komponenten zueinander in Beziehung gemaB der Gleichung , 



a > 

20 dann erhalt man einen Wert a, der ein MaB fiir die Nichtlinearitat ist. Zwar ware auch die 
Amplitute a 2 fur sich allein ein MaB fur die Nichtlinearitat, doch hangt a 2 auBer von der 
Nichtlinearitat auch noch von der Schwachung der durch die Anregung erzeugten Welle 
im Objekt ab. Der Parameter a ist davon weitgehend unabhangig. Er korrespondiert mit 
dem Wert a der Nichtlinearitat gemaB Gleichung 3. Diese Korrespondenz ist iiber die 

25 exakte Verwendung des Verzerrungstensors sowie der Freien Energie unter 

Berticksichtigung heherer Terme gegeben. Der Zusammenhang zwischen a und a ist 
streng monoton steigend. Deswegen kann analog zu Fig. 1 mittels Eo und a eine 
gleichwertige diagnostische Aussage gemacht werden. 
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Nachdem fur alle Voxel je ein Werte a ermittelt worden ist, lassen sich im Schritt 204 
Bilder erstellen, die die Nichtlinearitat der Elastizitat im Untersuchungsbereich 
darstellen. 

5 Da die Nichtlinearitat fur sich allein in der Regel noch keine Aussagen ttber die 
Gutartigkeit oder Bosartigkeit einer Lasion zulasst, wie in Verbindung mit Fig. 1 
erlautert, wird im Schritt 205 noch ein Bild erzeugt, das die raumliche VerteUung des 
Elastizitatsmoduls im Untersuchungsbereich darstellt. Die hierfiir erforderlichen 
Messwerte sind nach der Akquisition entsprechend Fig. 3 bereits vorhanden. Die 
10 Berechnung des Elastizitatsmoduls aus diesen Werten ist im Einzelnen in der eingangs 
erwahnten US-Appln. 09/743659 (PHD 99-152) dargestellt, worauf an dieser Stelle 
Bezug genommen wird. 

Der Benutzer kann dann einzelne, z. B. mittels eines Cursors markierte Bereiche in 
1 5 einem der Bilder herausgreifen und sich dafur die Werte von E und a - z. B. in einem 
Diagramm analog zu Fig. 1 - ausgeben lassen (Schritt 206). Danach endet das Verfahren 
207. 

Wie bereits ausgefiihrt, kann auch mit den Mitteln der Ultraschall-Elastographie die 
20 Ausbreitung der durch die mechanischen Schwingungen angeregten Wellen in dem 

Objekt verfolgt werden, und zwar wesentlich schneller als in der MR-Elastographie. Die 
aus den nichtlinearen Eigenschaften des Gewebes resultierenden nichtlinearen 
Verzerrung im zeitlichen Verlauf konnen dann in analoger Weise ermittelt werden, wie 
vorstehend erlautert. Neben der Mama konnen auch andere Bereiche untersucht werden, 
25 beispielsweise die Prostata. 
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PATENT ANSPRUCHE 



1. Verfahren zur Untersuchung eines Objektes mit den Schritten: 

a) Anregung des Objektes mit mechanischen Schwingungen, 

b) Messung des zeitlichen Verlaufs der Auslenkung von Voxeln des Objektes, die den 
von den Schwingungen im Objekt verursachten Weilen ausgesetzt sind, aus ihrer 

5 Ruhelage, 

c) Bestimmung der nichtlinearen Verzerrungen aus dem zeitlichen Verlaufs der . 
Auslenkung, 

d) Auswertung der nichtlinearen Verzerrungen. 

10 2. Verfahren nach Anspruch 1, wobei die Auswertung die Erzeugung eines die GroBe 
der nichtlinearen Verzerrungen darstellenden Bildes umfaBt 

3. Verfahren nach Anspruch 2, mit den Schritten: 

a) Anregung des Objektes mit zeitlich periodisch; vorzugsweise sinusfdrmig, 
15 verlaufenden mechanischen Schwingungen, 

b) Anregung der Kernmagnetisierung in dem Objekt in Verbindung mit einem zu den 
mechanischen Schwingungen synchronen magnetischen Gradientenfeld (Gl, G2) 
und Empfangen der entstehenden MR-Signale in dem Objekt zur Erzeugung eines 
MR-Phasen-Bildes wobei der zeitliche Verlauf des Gradientenfeldes so gewahlt ist, 

20 daB die MR-Signale sowohl von Auslenkungen mit der Grundfrequenz als auch 

von Auslenkungen mit wenigstens einer hdheren Harmonischen der Grundfrequenz 
bestimmt sind, 
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c) mehrfache Wiederholung des Schrittes b) unter Veranderung der Richtung des 
Gradienten des Gradientenfeldes und/oder der Phasendifferenz zwischen den 
mechanischen Schwingungen und dem Gradientenfeld zur Erzeugung weiterer 
Phasenbilder, 

d) . Bestimmung der Amplitude der Auslenkung der Kernspins mit der Grundfrequenz 

anhand der MR-Phasen-Bilder, 

e) Bestimmung der Amplitude der Auslenkung der Kernspins mit der hoheren 
Harmonischen der Grundfrequenz anhand der MR-Phasen-Bilder, 

f) Erzeugung eines vom Verhaltnis der Amplituden abhangigen Bildes. 



4. Verfahren nach Anspruch 3, mit den Schritten: 

wobei das magnetischen Gradientenfeld (Gl , G2) zwei einzelne sinusforrnige 
Schwingungen mit der halben Periodendauer der mechanischen Schwingungen umfaBt, 
deren zeitlicher Abstand voneinander ein Viertel der Periodendauer der mechanischen 
15 Schwingungen betragt. 

5. Verfahren nach Anspruch 4, wobei zwischen den beiden sinusfbrrnige Schwingungen 
ein refokussierender Hochfrequenzimpuls auf den Untersuchungsbereich einwirkt. 

20 6. Verfahren nach Anspruch 3 mit den Schritten: 

a) Zusatzliche Berechnung des Elastizitat im Untersuchungsbereich und 

b) Auswertung der fur die gleichen Voxel ermittelten Werte von Verzerrungen und 
Elastizitat. 

25 7. Anordnung zur Durchfuhrung des Verfahrens nach Anspruch 1 mit einem MR-Gerat, 
einem mechanischen Schwingungserzeuger, einer Auswerteeinheit, einem den zeitlichen 
Verlauf von magnetischen Gradientenfeldern bestimmenden Generator und einer das 
MR-Gerat, den Generator, den Schwingungserzeuger und die Auswerteeinheit 
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steuernden Steuereinheit, die so programmiert ist, daB folgende Schritte ausgeftthrt 
werden: 

a) Anregung des Objektes mit zeitlich periodisch, vorzugsweise sinusformig, 
5 verlaufenden mechanischen Schwingungen, 

b) Anregung der Kemmagnetisierung in dem Objekt in Verbindung mh einem zu den 
mechanischen Schwingungen synchronen magnetischen Gradientenfeld (Gl, G2) 
und Empfangen der entstehenden MR-Signale in dem Objekt zur Erzeugung eines 
MR-Phasen-Bildes wobei der zeitliche Verlauf des Gradientenfeldes so gewahlt ist, 

10 daB die MR-Signale sowohl von Auslenkungen mit der Grundfrequenz als auch 

von Auslenkungen mit wenigstens einer hoheren Harmonischen der Grundfrequenz 
bestimmt sind, 

c) mehrfache Wiederholung des Schrittes b) unter Veranderung der Richtung des 
Gradienten des Gradientenfeldes und/oder der Phasendifferenz zwischen den 

1 5 mechanischen Schwingungen und dem Gradientenfeld zur Erzeugung weiterer 

Phasenbilder, 

d) Bestimmung der Amplitude der Auslenkung der Kernspins mit der Grundfrequenz 
anhand der MR-Phasen-Bilder, 

e) Bestimmung der Amplitude der Auslenkung der Kernspins mit der hoheren 
20 Harmonischen der Grundfrequenz anhand der MR-Phasen-Bilder, 

f) Erzeugung eines vom Verhaltnis der Amplituden abhangigen Bildes 

8. Computerprogramm fur eine auf ein MR-Gerat, einen Schwingungserzeuger und eine 
Auswerteeinheit einwirkende Steuereinheit zur Durchfuhrung des Verfahrens nach 
25 Anspruch 1 gemaB folgendem Ablauf: 

a) Anregung des Objektes mit zeitlich periodisch, vorzugsweise sinusformig, 
verlaufenden mechanischen Schwingungen, 

b) Anregung der Kemmagnetisierung in dem Objekt in Verbindung mit einem zu den 
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mechanischen Schwingungen synchronen magnetischen Gradientenfeld (Gl, G2) 
und Empfangen der entstehenden MR-Signale in dem Objekt zur Erzeugung eines 
MR-Phasen-Bildes wobei der zeitliche Verlauf des Gradientenfeldes so gewahlt ist, 
dafi die MR-Signale sowohl von Auslenkungen mit der Grundfrequenz als auch 
von Auslenkungen mit wenigstens einer haheren Harmonischen der Grundfrequenz 
bestimmt sind, 

mehrfache Wiederholung des Schrittes b) unter Veranderung der Richtung des 
Gradienten des Gradientenfeldes und/oder der Phasendifferenz zwischen den 
mechanischen Schwingungen und dem Gradientenfeld zur Erzeugung weiterer 
Phasenbilder, 

Bestimmung der Amplitude der Auslenkung der Kernspins mit der Grundfrequenz 
anhand der MR-Phasen-Bilder, 

Bestimmung der Amplitude der Auslenkung der Kernspins mit der hOheren 
Harmonischen der Grundfrequenz anhand der MR-Phasen-Bilder, 
Erzeugung eines vom Verhaltnis der Amplituden abhangigen Bildes 
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ZUS AMMENF AS SUNG 

Verfahren zur elastographischen Untersuchung eines Objektes 

Die Erfindung betrifft ein Elastographie-Verfahren, dessen diagnostische Aussagekraft 
dadurch verbessert wird, dass der zeitliche Verlauf der Auslenkung, mit dem ein durch 
mechanische Schwingungen angeregter Bereich reagiert, analysiert wird und die 
nichtUnearen Verzerrungen registriert werden. Diese Verzerrungen sind ein MaB fur die 
nichtlinearen elastischen Eigenschaften des Bereiches und stellen eine diagnostisch 
relevante Information dar. 



10 Fig. 4 
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